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SUMMARY 


The flutter characteristics of two untapered wings, one having 15° of sweep- 
back and an aspect ratio of 5*34 and the other having 30 ° of sweepback and an 
aspect ratio of 4.1 6 , have been investigated analytically for Mach numbers up to 
3.0. The calculations employed piston theory or quasi-steady second-order theory 
and included variations in the number and type of vibration modes employed (i.e., 
measured natural modes or uncoupled beam-type modes) and in the treatment of 
finite wing thickness and of loading in the tip region. Results of the analysis 
and comparison with experiments indicated that the accuracy of flutter speeds 
predicted by piston theory and by quasi-steady second-order theory varied consid- 
erably, depending on several factors such as the number and type of vibration 
modes employed (coupled or uncoupled) and on the treatment of the aerodynamic 
loading near the streamwise wing tips. Best results for the two wings were not 
obtained by the same combination of these factors. 


INTRODUCTION 


In numerous flutter analyses for high supersonic Mach numbers, piston theory 
(refs. 1 and 2) or quasi-steady second-order theory (refs. 3 and 4) has been 
employed for the calculation of the required oscillatory aerodynamic forces. The 
pressure expressions resulting from these two theories are simple to evaluate and 
are well suited for flutter-trend studies and theoretical comparisons at very high 
Mach numbers. Many of the applications, however, have been for theoretical trend 
studies of two-dimensional wing sections (for example, refs. 5 to 7) or for com- 
parisons with other flutter analyses (for example, refs. 8 and 9)- On the other 
hand, comparisons with experimental flutter data for swept wings of finite span 
at moderate supersonic Mach numbers do not appear to be extensive, although some 
results are shown in references 4 and 10. In the flutter analysis of finite wing 
planforms, the effects and interactions of several pertinent factors require fur- 
ther examination. These factors include Mach number, leading-edge sweep, finite 
airfoil thickness, loading near wing tips, and number and type of vibration modes. 



The influences of these factors on flutter characteristics calculated by 
second-order piston theory and by quasi- steady second-order theory are investi- 
gated herein for two homogeneous untapered swept wings at Mach numbers from 1.65 
to 5.00. The results are compared with the experimental flutter data of refer- 
ence 11. Flutter calculations were made both with and without consideration of 
finite wing thickness. In these calculations three different methods were used 
for evaluating the aerodynamic loading within the triangular region affected by 
the streamwise wing tips. The calculations were performed with two or three natu- 
ral (coupled) vibration modes or with two, three, or four uncoupled vibration 
modes. 
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reduced frequency, —y— 

Mach number 

total mass of wing panel 

lifting pressure 

time 

flutter speed 

volume of air within a cylinder having the streamwise chord as base 
diameter and the panel span as height 

streamwise distance measured rearward from wing leading edge, fraction 
of local chord 

spanwise distance measured from wing root, fraction of panel span 

magnitude of airfoil ordinate measured normal to chord plane, fraction 
of local chord 

amplitude of vertical displacement at point (x,y) on wing panel 
modal deflection for k^ coupled vibration mode 
modal deflection for i^ uncoupled bending vibration mode 
modal deflection for uncoupled torsional vibration mode 


p = \Jbi 2 - 1 

y ratio of specific heats 
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pv 


sweep angle 

mass ratio for wing panel, 

air density 
circular frequency of vibration at flutter 
circular frequency of coupled vibration mode 

circular frequency of i ^ 1 uncoupled bending vibration mode 
circular frequency of uncoupled torsional vibration mode 


WINGS 


General Description 

The flutter calculations of this report were made for two of the untapered, 
swept wings for which experimental flutter data were given in reference 11. One 
wing was swept back 15 ° and had an aspect ratio of 5 - 3 ^; the other wing was swept 
back 30° and had an aspect ratio of k.l6. (See fig. 1.) These wings were cut 
from 0. 04l-inch- thick magnesium sheet, and had chords of 2 inches measured perpen- 
dicular to the leading edge. As shown in figure 1, the leading and trailing edges 
of both wings were beveled l/4 inch to form a 2 . 05 ~percent-thick symmetrical hex- 
agonal airfoil section perpendicular to the leading edge. For these homogeneous 
untapered wings, both the elastic axis and the local centers of gravity were taken 
to be at midchord. 

In the investigation of reference 11, two magnesium models and two aluminum 
models of each planform were flutter tested. For convenient reference, the models 
that were tested at supersonic speeds are designated as A, B, and C. Flutter 
calculations are presented herein only for 15° wing models A and B, and for 30° 
wing models A and B (see table I), although experimental flutter data for 15° 
and 30° wing models C are included to indicate the experimental trends. 

These wings seemed particularly appropriate for this investigation for 
several reasons: 

(1) Because of the sharp leading edges and relatively small sweep angles, 
piston theory and quasi- steady second-order theory should be reasonably accurate 
at fairly low supersonic Mach numbers. 

(2) The large streamwise tips should have a sizable effect on the aero- 
dynamic loading. 

( 3 ) Some of the measured natural vibration modes for these wings (ref. 12) 
indicate a high degree of bending- torsion coupling so that coupled-mode and 
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uncoupled-mode representations of these wings should provide quite different 
inputs to the flutter analysis. 

(4) Experimental flutter data (ref. 11 ) were available for comparison. 


Mode Shapes and Frequencies 

The first three natural (coupled) mode shapes for both wings were measured 
by a bouncing sand technique. The resulting deflections were presented in 
reference 12 and are repeated for convenience herein in figures 2(a) and 3( a )* 

The corresponding measured coupled-mode frequencies for the models of this inves- 
tigation were given with the flutter test data in reference 11, and are summarized 
in table I of the present report. The node lines for these wings (fig. l) indi- 
cate that the second and third natural modes involve a high degree of bending- 
torsion coupling, particularly for the 30 ° wing. 

Since the wings employed herein were homogeneous and untapered, the uncoupled 
mode shapes and frequencies for both planforms were considered to be the same as 
those for a uniform cantilever beam. The use of uncoupled modes implies that 
wing sections normal to the elastic axis (midchord for the present two wings) are 
considered to oscillate without distortion. For swept wings, these sections do 
not lie in the f ree-stream direction, so that some camber of streamwise sections 
appears (figs. 2(b) and 3(t>)) as it does in the measured coupled modes. The 
first three bending-mode and first torsion-mode frequencies for 15° wing models A 
and C and the first two bending-mode and first two torsion-mode frequencies for 
30° wing models A and C were calculated from the simple formulas for a uniform 
cantilever beam. To allow for slight differences in stiffness between ostensibly 
identical models A and B, the uncoupled-mode frequencies for model B were obtained 
from those for model A by multiplying by the ratio of the measured first-mode 
frequency of model B to the corresponding frequency of model A; that is, for 
models B of both planforms, 


(^u)b = psjjj x K^a 

and. 

M 

(“ a > 4 b = x (^J)a 

The uncoupled-mode frequencies for all models are presented in table I. For each 
model, the frequency of the second bending mode is close to that of the first 
torsion mode, and both are near the measured second- and third-mode natural 
(coupled) frequencies. 
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FLUTTER CALCULATIONS 


The aerodynamic forces for the present flutter calculations were based on 
the expressions for lifting pressure given in reference 4 for second-order piston 
theory and for quasi-steady second-order theory. The quasi-steady second-order 
theory employed herein and in reference 4 is based on the steady-flow second-order 
theory of reference 3 and should not be confused with the unsteady second-order 
theory of reference 13. The quasi-steady theory was used here because the meas- 
ured reduced frequencies at flutter for the present wings were small (table i) 
and because the calculations could be made with the same computing program used 
for the piston-theory calculations. Accordingly, for harmonic motion the expres- 
sions for lifting pressure may be written as 


&p(x,y, t) _ 
pV 2 /2 

for piston theory and as 
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for quasi-steady second-order theory. These lifting- pressure expressions are 
similar in form, with second-order theory differing from piston theory only in two 
coefficients that are functions only of Mach number and the ratio of specific 
heats. Both theories permit the inclusion of airfoil shape and coupled vibration 
modes with camber. 

In the quasi-steady theory (ref. 4) the pressure expression in equation (2) 
is used to represent the pressure distribution on an oscillating wing as though 
it were composed of a succession of steady-state distributions each associated 
with an instantaneous angle of attack. This procedure is reasonable for high 
Mach numbers because the reduced frequency at flutter generally decreases as Mach 
number increases. Thus, at high Mach numbers the unsteady aspects of the flow 
have reduced significance. 

Most of the flutter calculations for the wings of this report included the 
effect of finite wing thickness and employed the first three coupled vibration 
modes or three uncoupled modes (first and second bending and first torsion). 

Some additional calculations for both wings neglected finite wing thickness, and 
some employed only the first two coupled modes or two uncoupled modes (first 
bending and first torsion) for comparison purposes-*-. (See table II. ) The fourth 

-*-In view of the closeness of the second and third natural frequencies and of 
the first torsion and second bending frequencies for these two wings (table I) 
and in view of the high degree of bending-torsion coupling shown in both the 
second and third natural modes (ref. 12), the use of at least three vibration 
modes would ordinarily be indicated for flutter analysis of these wings. The 
two-mode calculations included herein are therefore considered to be of analytical 
interest only. 
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uncoupled mode (third bending for the 15 ° swept wing, second torsion for the 

5 wept wing) was included in some uncoupled-mode calculations in order to indi- 
cate the degree of convergence of the present modal representations. The coupled 
modes were not assumed orthogonal, and the off-diagonal generalized masses, though 
small, were included in the flutter determinant. For the uncoupled-mode calcula- 
tions, however, the off-diagonal generalized masses were identically zero. . Since 
the pertinent structural damping values are not known for the models used in this 
investigation, and since the damping coefficients for cantilevered solid metal 
wings of the present type are usually very small, all calculated flutter points 
herein are taken to be points for which the required structural damping vanishes. 
All integrations required for the evaluation of the generalized aerodynamic forces 
and of the generalized masses were performed by use of an 11-point Simpson's rule. 

Piston theory and quasi-steady second-order theory as expressed in refer- 
ence 4 take no account of streamwise wing tips. However, for some of the calcula- 
tions of this report and of reference l4, an approximate tip- load modification was 
made on the basis of steady-flow linear theory. The streamwise tip, of course, 
influences loading only within the triangular region bounded by the tip, the 
trailing edge, and the Mach line from the leading-edge tip; and the area of this 
triangle decreases as Mach number increases. Furthermore, the reduced frequency 
at flutter generally decreases as Mach number increases so that a steady-flow type 
of tip- load modification should be more accurate at the higher Mach numbers. The 
tip- load modification as applied herein and in reference l4 consists of multi- 
plying the piston- theory or second-order- theory loading at each point within the 
tip triangle by the ratio of steady-state load with streamwise tip to steady-state 
load without streamwise tip. This ratio, of course, varies with the location of 
the point within the tip triangle and with Mach number. Thus for piston theory, 
for example, the corrected lifting pressure at a point (x,y) on the wing is given 

"by 


Ap(x,y,M) piston theory corrected 


Ap(x,y,M) plston theory uncorrected X Ap(x,y,M) 


Ap(x,y,M) linear p heor y with tip 


(3) 


linear theory without tip 


Such calculations for the wings of this report should give an indication of 
the merit of this modification because these wings were untapered and of moderate 
aspect ratio so that the tip triangle included a sizable portion of the wing area 
at the lower Mach numbers. The sensitivity of the calculations to this modifica- 
tion was examined by computing flutter characteristics for the 15 wing with no 
tip- load modification, with the above-mentioned tip- load modification, and with 
zero load assumed within the tip triangle. (See table II.) The first and last 
of these procedures represent limiting cases. . For the 30° wing, either no tip- 
load modification or the calculated modification was used. 
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PRESENTATION OF RESULTS 


A summary of pertinent conditions for all the calculations and an index to 
the results are provided in table II. In figures 4 to 15 the supersonic flutter 
speeds and frequencies calculated for the 15 ° and 50 ° wings by piston theory and 
by quasi-steady second-order theory are compared with the flutter experiments of 
reference 11. Figures 4 to 7 and 10 to 15 show the effects on the calculated 
flutter characteristics of the number and type of vibration modes employed and of 
the inclusion of finite wing thickness. Figures 8 and 9> together with fig- 
ures 10 , 11 , 14, and 15 , present the results obtained with the three different 
representations of loading in the tip region. 


DISCUSSION OF RESULTS 


In the calculations of all the flutter-speed and flutter- frequency curves of 
figures 4 to 15 , the wing properties and flow density employed were those associ- 
ated with the experimental flutter points shown. Since each experimental point 
(from ref. ll) represents a different model and a different density, the calcula- 
tions for each model and density were limited to a range of Mach number which 
bracketed the experimental value. Specifically, the Mach numbers covered are 
2.50 to 5«00 for model A and I .65 to 2.50 for model B. Because of the discontin- 
uous change of mass ratio and modal frequencies, the calculated flutter character- 
istics show discontinuities at M = 2.50. 

Figures 4 to 15 indicate that for the present two wings the closest agree- 
ment between theory and experiment was not attained with the same combination of 
the variable factors (i.e., piston theory or quasi-steady second-order theory, 
number and type of vibration modes, wing thickness, and tip- load modification). 
Furthermore, the effects of varying these factors were frequently different for 
the two wings. The subsequent sections of this report discuss the results 
obtained by changing each of these factors individually. 


Comparison of Results From Piston Theory and From 
Quasi-Steady Second-Order Theory 

The flutter speeds calculated for these wings from quasi-steady second-order 
theory were lower than those calculated from piston theory. In addition to being 
lower, the flutter-speed coefficients calculated by second-order theory generally 
followed the experimental trend of flutter-speed coefficient with Mach number more 
closely, particularly at the lower supersonic Mach numbers . Furthermore, there 
was usually little difference between the flutter frequencies obtained from these 
two theories. Since the lifting- pressure expressions associated with piston the- 
ory and with quasi-steady second-order theory (ref. 4) approach each other as Mach 
number approaches infinity, the flutter speeds indicated by the two theories also 
approach each other as Mach number increases. In view of these statements, the 
following discussion is concerned with the results of the second-order- theory 
calculations unless otherwise specified. 
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Results for 15° Wing 


Humber and type of vibration modes .- For the 15° wing, results of the flutter 
calculations which include two, three, or four uncoupled vibration modes, finite 
wing thickness, and no modification of aerodynamic loading in the tip region are 
shown in figures 4(a) and 5(a). For all Mach numbers the calculated flutter- 
gpgg(j coefficients decrease as the number of uncoupled modes is increased. How- 
ever, the difference between three-mode and four-mode calculations is insignifi- 
cant except at the lowest Mach numbers where the four-mode piston- theory calcula- 
tion produced two flutter solutions which indicated that flutter could occur in 
either of two modes. Only the solution yielding the lower flutter speed is shown 
in figures 4(a) and 5(a). The corresponding second-order- theory calculations 
indicated the possibility of a similar double solution at Mach numbers lower than 
those calculated. 

Figures 4(b) and 5(b) present flutter characteristics for the 15° wing cal- 
culated with two or three coupled vibration modes, finite wing thickness, and no 
tip- load modification. In contrast to the results with uncoupled modes, the addi- 
tion of the third coupled mode increases the calculated flutter speed at all Mach 
numbers. Since higher coupled modes are not available from the measurements of 
reference 12, it is not known whether their inclusion in the flutter calculations 
would lead to further changes in the results. 

The flutter speeds calculated with three coupled modes (fig. 4(b)) are lower 
than those calculated with three uncoupled modes (fig. 4(a)). (See also 
fig. 8(a).) For all Mach numbers, however, the difference is less than 10 per- 
cent. For these two calculations, the agreement with experiment is satisfactory 
at Mach number 3-0, but becomes poorer as Mach number decreases. 

The calculated flutter frequency ratios (figs. 5(a) and (b)) are in satisfac- 
tory agreement with experiment and generally show little change as the number of 
vibration modes is increased. 

Calculations similar to those represented in figures 4 and 5 but without 
finite wing thickness yield the flutter speeds and frequencies shown in figures 6 
and 7. Inclusion of the third uncoupled mode (fig. 6(a)) depresses the calculated 
flutter speed, whereas inclusion of the third coupled mode (fig. 6(b)) raises the 
flutter speed. This result indicates changes that are in the same direction as 
when thickness is included. The magnitudes of the changes in figure 6 (thickness 
neglected) are much greater than those in figure 4 (thickness included). For 
example, the two-uncoupled-mode calculation gave no solution when thickness was 
omitted. 

Wing thickness.- Flutter calculations employing two uncoupled vibration _ modes 
yield finite flutter speeds when finite wing thickness is included (fig. 4(a)), 
but not when thickness is neglected. Hence, in this case the effect of thickness 
is large and destabilizing. This behavior is in qualitative agreement with the 
results of two-dimensional piston- theory calculations for low ratios of bending 
frequency to torsion frequency (refs. 5 and. 7 to 9)- I n reference 15 flutter 
experiments and representative-section piston- theory calculations for a wing of 
rectangular planform also showed a destabilizing effect of thickness. In contrast 
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to these results, the omission of finite wing thickness in three-uncoupled-mode 
calculations is destabilizing and yields conservative flutter speeds* (Compare 
figs. 4(a) and 6(a).) The binary flutter trends for two-dimensional wings in fig- 
ure 2 of reference 9 and in figure 13 of reference 7 also show that the inclusion 
of finite thickness can be stabilizing in some cases. The appreciable difference 
in levels of flutter-speed coefficient between models A and B (see fig. 6(a) at 
M = 2.50) is caused entirely by the difference in flow density shown in table I; 
thus, for these conditions the calculated flutter dynamic pressure is not invari- 
ant with changes in density. 

The inclusion of finite thickness in flutter calculations with two coupled 
vibration modes (compare fig. 4(b) and 6(b)) has a destabilizing effect that is 
small in comparison with that indicated by the calculations with two uncoupled 
modes even though the modal-frequency ratios are small in both cases. On the 
other hand, in the calculations with three coupled modes, inclusion of thickness 
has large depressing effects on both flutter speed and frequency. (See 
figs. 4(b), 5(b), 6(b), and 7(h).) Similar calculations with three coupled modes 
for an untapered 45° swept wing (ref. 4) showed little effect of wing thickness. 
Since the airfoil (normal to the leading edge) for this 45° swept wing was the 
same as that for the present 15° and 30° wings, the thickness ratio and contour 
slopes in the streamwise direction were lower for the wing of reference 4. Fur- 
ther indications of the influence of wing thickness and other parameters on 
representative-section piston-theory flutter trends are presented, for example, 
in references 5 ho 8. 

It is evident from figures 6(b) and 7(h) that the calculations with three 
coupled modes and no thickness are far more sensitive to changes in flow density 
and modal- frequency ratios than are any of the previously mentioned calculations. 
At M = 2.50, the only quantities differing in the coupled-mode calculations for 
models A and B are the flow density and the ratio of second-mode to third-mode 
frequency. (See table 1(a).) Some further discussions of the* effects of flow 
density are included in references l4 and l6. 

Load representation in the tip region .- Progressive reduction of the aero- 
dynamic loading in the triangular region influenced by the tip from the no-tip- 
modification level to the calculated tip-modification level and then to the zero- 
load condition is shown in figures 8 and 9 ho cause progressive reduchion in the 
calculated flutter- speed coefficients and flutter frequency ratios. Furthermore, 
the use of the calculated tip- load modification (calculated from steady-flow the- 
ory as previously described) reduces the calculated flutter speed more at lower 
Mach numbers than at higher Mach numbers. (Compare figs. 8(a) and 8(b).) Hence, 
when the calculated tip- load modification is employed, the slope of the curve of 
flutter-speed coefficient against Mach number (fig. 8(b)) increases and corre- 
sponds more closely to the experimental trend. In comparison, the limiting case 
of zero load in the tip region (fig. 8(c)) gives still lower flutter speeds, but 
the trends appear increasingly conservative as Mach number increases. In all 
cases, the lowest flutter speeds and frequencies are obtained from quasi-steady 
second-order theory with coupled vibration modes. In contrast, in similar calcu- 
lations for the highly tapered 45° swept wing of reference l4, the lowest flutter 
speeds were obtained from the second-order theory with uncoupled vibration modes. 
However, use of the calculated tip- load modification in reference 14 lowered the 
calculated flutter speed as in the present case. 
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Best results .- For the 15° wing, the best overall predictions of flutter- 
speed trends and ma gn itudes were obtained from quasi-steady second-order- theory 
calculations in which finite wing thickness was included and in which coupled 
vibration modes and the calculated tip- load modification were employed 
(fig. 8(b)). Even for this calculation the predicted flutter speed at M = 2.00 
is unconservative by about 10 percent. For comparison, the rule of thumb given 
in reference 5 for the applicability of piston theory, namely that M cos A ^ 2.5? 
is satisfied for the 15*^ wing only for M ^ 2.6. Above this Mach number the 
agreement between the present best results and experiment appears to be very good. 


Results for 50° Wing 

Number and type of vibration modes . - For the 50° wing, results of the flutter 
calculations which include two, three, or four uncoupled vibration modes, finite 
wing thickness, and no modification of aerodynamic loading in the tip region are 
shown in figures 10 and 11. The flutter-speed coefficients and flutter frequency 
ratios calculated with two modes are excessively high and no flutter solutions 
occurred for M ^ 2.40. In contrast, inclusion of the second bending mode yields 
flutter speeds that are slightly conservative and in good agreement with experi- 
mental values, although the calculated flutter frequencies remain high. Thus, the 
second bending mode has a much more significant effect on the flutter results for 
the 50° wing than for the 15° wing. (Compare figs. 4(a) and 10 and figs. 5(a) 
and 11.) As in the case of the 15° wing, the inclusion of the fourth uncoupled 
mode, which for the 50 ° wing is the second torsion mode, gives flutter speeds and 
frequencies which are essentially the same as the three-mode values. 

When finite wing thickness is neglected (figs. 12 and 15), the results for 
the 50° wing are qualitatively the same as for the 15 ° wing (figs. 6(a) and 7(a))- 
That is, no flutter solution occurred when two modes were used, whereas the flut- 
ter speeds obtained with three modes were excessively conservative. 

Similar calculations employing two or three coupled vibration modes, with or 
without finite wing thickness and without the tip- load modification, gave no flut- 
ter solution. In contrast, similar calculations for the 15° wing of the present 
report and for the untapered 45° wing of reference 4 yielded finite flutter speeds 
and frequencies. The calculated flutter speeds for the wing of reference 4, how- 
ever, were unconservative, and the calculated flutter frequencies were nearly 
twice the experimental values. 

Wing thickness . - Figures 10 and 12 indicate that when uncoupled vibration 
modes are employed, changes in the calculated flutter speed caused by the inclu- 
sion of wing thickness are in the same direction as those shown previously for 
the 15° wing. In the two-mode calculations, flutter solutions occurred only at 
the higher Mach numbers and only when thickness was included. Therefore, the 
effect of thickness was large and destabilizing, at least at these higher Mach 
numbers. The inclusion of thickness in the three-mode calculations again was 
stabilizing, particularly at the higher Mach numbers. 

Load representation in the tip region . - Comparisons of figures 10 and 14 and 
figures 11 and 15 show the effects of employing the calculated modification to the 
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aerodynamic loading in the tip region for the calculations with three vibration 
modes and with finite thickness. For the 30° wing, no calculations were made for 
the limiting condition of zero load in the tip region. At the higher Mach num- 
bers, the uncoupled-mode calculations for the 30 ° wing show that use of the calcu- 
lated tip-load modification causes a moderate reduction in flutter speed as it did 
in the case of the 15° wing. At Mach numbers near 2.00, however, the tip- load 
modification has little effect, and at the lowest Mach numbers it even increases 
the calculated flutter speed (figs. 10 and Ik). These varying effects of the tip- 
load modification with Mach number may be related to two predominant influences 
which the modification has in the flutter analyses. First, application of the 
calculated tip- load modification reduces the magnitude of the aerodynamic loading 
in the tip triangle; this reduction would tend to increase the calculated flutter 
speed. Second, use of the tip- load modification shifts the center of loading 
forward in the tip region; this shift generally tends to decrease the calculated 
flutter speed. Apparently, for the 30° wing, the former effect dominates at the 
lower Mach numbers, and the latter effect dominates at the higher Mach numbers. 

In calculations employing uncoupled vibration modes, use of the tip- load modifica- 
tion caused negligible change in the calculated flutter frequencies which remain 
considerably higher than experimental values (fig. 15)- 

Although no coupled-mode flutter solution was obtained without the tip- load 
modification, use of the modification in coupled-mode flutter calculations yielded 
finite flutter speeds; however, these results are unconservative and are not in 
good agreement with experiment (fig. lk) . The corresponding calculated flutter 
frequencies are in fair agreement with experiment (fig. 15)- These results indi- 
cate that coupled-mode calculations for the 30 ° wing are highly sensitive to the 
treatment of the loading in the tip region. 

Best results .- For the 30° wing, the combination of factors yielding the best 
calculated flutter characteristics is not as clearly defined as for the 15 ° wing. 
Piston-theory and quasi-steady second-order- theory flutter calculations employing 
three uncoupled vibration modes and finite wing thickness, with or without the 
calculated tip- load modification, give flutter speeds that are close together and 
moderately conservative at all but the lowest Mach numbers. Since these flutter 
speeds are conservative and since the values given by piston theory are slightly 
higher than those given by the quasi- steady second-order theory, the piston-theory 
values are slightly closer to the experimental points. On the other hand, more 
accurate values of flutter frequency are obtained from the second-order theory 
with three coupled modes, finite thickness, and the calculated tip- load modifica- 
tion (fig. l4), although the corresponding flutter speed is about 50 percent high 
at M = 2.00. The rule of thumb for the applicability of piston theory (ref. 5)> 
namely that M cos A ^ 2.5> is satisfied for the 30° wing only for M ^ 2.9- In 
view of this criterion and the previous discussion, the most satisfactory calcula- 
tions for the 30 ° wing are considered to be obtained from either piston theory or 
quasi-steady second-order theory with three uncoupled vibration modes, finite 
wing thickness, and no tip- load modification. 
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CONCLUDING REMARKS 


The flutter characteristics of two untapered swept wings have been investi- 
gated analytically for Mach numbers up to 3.0 by piston theory and by quasi-steady 
second-order theory. The calculations included variations in the number and type 
of vibration modes employed (i.e., measured natural modes or uncoupled beam-type 
modes) and in the treatment of finite wing thickness and of loading in the tip 
region. 

For the present two wings, the closest overall agreement between theory and 
experiment was not attained with the same combination of the previously mentioned 
variable factors, and the effects of varying these factors were frequently differ- 
ent for the two wings. For example, the best results for the 15° wing were 
obtained from the quasi-steady second-order theory with coupled vibration modes, 
finite wing thickness, and a modification of the loading near streamwise wing 
tips based on three-dimensional steady-flow theory. On the other hand, for the 
30 ° wing, flutter characteristics given by piston theory and by quasi- steady 
second-order theory were generally close together, and best results were obtained 
with uncoupled vibration modes, finite wing thickness, and no modification of 
loading in the tip region. Furthermore, the effect of including either finite 
wing thickness or the calculated tip- load modification can be large or small and 
either stabilizing or destabilizing depending, for example, on number and type of 
vibration modes, modal- frequency ratios, planform, and Mach number. However, the 
quasi- steady second-order theory generally yielded calculated flutter speeds which 
were lower and in better agreement with experimental trends than those given by 
piston theory, particularly at the lower Mach numbers. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 8, 1963* 
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TABLE I.- MODAL FREQUENCIES AND FLUTTER TEST CONDITIONS 
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1,319 1,721 267 1,673 1,533 
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Figure k.- Calculated flutter-speed coefficients for 15° wing shoving effects of number of vibration modes employed. 
All calculations include finite wing thickness and no modification of load in the tip triangle. 




No. modes Aerodynamic theory 
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(b) Coupled modes. 



o 
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Figure Calculated flutter-frequency ratios for 15° wing showing effects of number of inodes employed. All calculations 
include finite wing thickness and no modification of load in tip triangle. 



o 
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Figure 5-- Concluded. 



Flutter calculations 



Figure 6.- Calculated flutter-speed coefficients for 15° wing showing effects of number of vibration modes employed 
All calculations neglect finite wing thickness and include no modification of load in tip triangle. 














(a) No modification of load in tip triangle. 

Figure 8.- Calculated flutter-speed coefficients for 15° wing with three different load 
representations in the tip region. All calculations include finite wing thickness 
and three vibration modes. 
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(b) Calculated modification of load in tip triangle. 
Figure 3 . - Continued . 
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Figure 8.- Concluded. 
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Modes Aerodynamic theory 
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(a) No modification of load in tip triangle. 

Figure 9-- Calculated flutter- frequency ratios for 15° wing with three different load 
representations in the tip region. All calculations include finite wing thickness 
and three vibration modes. 
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Flutter calculations 







Modes Aerodynamic theory 
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(b) Calculated modification of load in tip triangle. 
Figure 9.- Continued. 
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(c) Zero load in tip triangle. 
Figure 9-- Concluded. 
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Figure 10.- Calculated flutter- speed coefficients for 30° wing shoving effects of number of 
uncoupled vibration modes employed. All calculations include finite wing thickness and 
no modification of load in tip triangle. 






Flutter calculations 



Figure 12.- Calculated flutter-speed coefficients for J0° wing with uncoupled vibration modes. These calculations 
neglect finite thickness and include no modification of load in tip triangle. 




Figure 13*- Calculated flutter- frequency ratios for 30° wing with uncoupled vibration modes. These calculations neglect 
finite wing thickness and include no modification of load in tip triangle. 




Flutter calculation# 



Figure lh.- Calculated flutter-speed coefficients for 30° wing. All calculations employ calculated modification of load 

in tip triangle, finite wing thickness, and three vibration modes. 
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Figure 15.- Calculated flutter- frequency ratios for 30° wing. All calculations employ calculated modification of load in 

tip triangle, finite wing thickness, and three vibration modes. 





